AEI-2009-062 
ULB-TH/09-23 



Supersymmetric quantum cosmological billiards 



Axel Kleinschmidt 
Physique Theorique et Mathematique & International Solvay Institutes, 
Universite Libre de Bruxelles, Boulevard du Triomphe, ULB-CP231, BE-1050 Bruxelles, Belgium 



o 
o 

oo 



cr 

i , 



(N 
> 
00 

o 

rn 
o 

On 
O 



X 



Michael Koehn and Hermann Nicolai 
Max-Planck-Institut fur Gravitationsphysik, Albert- Einstein- Institut, Am Muhlenberg 1, DE-14476 Golm, Germany 

D = 11 Supergravity near a space-like singularity admits a cosmological billiard description based 
on the hyperbolic Kac-Moody group Eio. The quantization of this system via the supersymmetry 
constraint is shown to lead to wavefunctions involving automorphic (Maass wave) forms under the 
modular group W + (E\o) = PSL<2(0) with Dirichlet boundary conditions on the billiard domain. 
A general inequality for the Laplace eigenvalues of these automorphic forms implies that the wave 
function of the universe is generically complex and always tends to zero when approaching the initial 
singularity. We discuss possible implications of this result for the question of singularity resolution 
in quantum cosmology and comment on the differences with other approaches. 
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One of the surprising results in the study of classical 
gravity very close to a generic space-like singularity due 
to Belinskii, Khalatnikov and Lifshitz (BKL) [l[ was the 
realization that, under rather general assumptions, the 
system becomes ultralocal in space and can be described 
by a sequence of Kasner regimes. In the strict limit to- 
wards the singularity, the Kasner behavior is interspersed 
with hard reflections of the logarithms of the spatial scale 
factors off infinite potential walls [2, 3]. This behavior has 
been termed 'cosmological billiards'; the geometry of the 
billiard table as well as the possible occurrence of chaotic 
oscillations near the singularity depend on the dimension 
and the matter content of the theory [J, [f| . In particu- 
lar, for D = 11 supergravity it was shown by Damour 
and Henneaux Q that the billiard domain is the funda- 
mental Weyl chamber C of the hyperbolic Kac-Moody 
group Eio, implying chaotic behavior. 

Resolving the cosmological singularity requires the 
transition to the quantum theory. In this paper we ad- 
dress this issue in the framework of cosmological billiards. 
More precisely, we will set up and solve the quantum 
constraints for D = 11 supergravity for the ten spa- 
tial scale factors and the fermionic degrees of freedom 
in compliance with the supersymmetry constraint. The 
approach followed here is thus a variant of 'minisupcr- 
space' quantization of gravity pioneered in 0, H, 0] and 
further developed in [ljl El El El 0, El • An essen- 
tial new ingredient of the present work is the arithmetic 
structure provided by E\o and its Weyl group, whose rel- 
evance in the context of Einstein gravity was pointed out 
and explored in [16J. In accordance with the M-theory 
proposal of [I?) , where a correspondence was established 
at the classical level between a (truncated) gradient ex- 
pansion of the D = 11 supergravity equations of motion 
and an expansion in heights of roots of a constrained 
'geodesic' Eio/K(Eio) coset space model, the cosmolog- 
ical billiards approximation corresponds to the restric- 
tion of the coset model to the Cartan subalgebra of E w . 
Our results therefore represent the first step towards the 



quantization of the full coset model. 

Starting with the bosonic variables, the diagonal met- 
ric considered for the cosmological billiard for a (d + 1)- 
dimensional space-time is of the form [l8[ (for d > 3) 
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leading to the kinetic term £kin = \n 1 Yl a b=i P a G a bl3 b 
in terms of the lapse n = Ng -1 (the spatial volume is 
g 1 / 2 = exp[— ~J2 a P a }) and the Lorentzian DeWitt met- 
ric [m 
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It will be essential that for d = 10 this metric coincides 
with the restriction of the Cartan-Killing metric of Eiq 
to its Cartan subalgebra. The spatial ultralocality of 
the BKL limit reduces the gravitational model to a clas- 
sical mechanics system of a relativistic billiard ball de- 
scribed by the [3 a variables moving on straight null lines 
in the Lorentzian space with metric G a b until hitting a 
billiard table wall. The straight line segments are Kasner 
regimes. The singularity is at t = +oo in the 'Zeno-likc' 
time co-ordinate t that is related to physical (proper) 
time T by t ~ — logT. There is one such system for each 
spatial point x, and these systems are all decoupled. 

The conjugate canonical bosonic variables of the bil- 
liard system are (3 a and 7r a — G a bf3 b ; the Hamiltonian 
then is H = ^n a G ah Hb with the inverse metric G ab , and 
we have set the lapse n = 1. Before quantisation, we 
perform the following change of variables by means of 
which the billiard motion can be projected onto the unit 
hyperboloid in /3-space 0| 



P = pW , LU GabOO 



1 . 



p 2 = -l3 a G ab l3 b , (3) 



where p is the 'radial' direction in the future light-cone 
and oj a = u> a (z) are expressible as functions of d — 1 
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coordinates z on the unit hyperboloid. The limit towards 
the singularity is p — > oo in these variables. The Whcclcr- 
DeWitt (WDW) operator on /3-space takes the form [2(3] 

n = G ab a a d b = - P l - d §- p (p d ^£) + ^' A lb , (4) 

where Alb is the Laplace-Beltrami operator on the 
(d — l)-dimensional unit hyperboloid. The WDW equa- 
tion therefore reads H&(p, z) — for the wavefunction 
$(p, z). As usual (see e.g. [15|]) one can adopt p as a time 
coordinate in the initially 'timeless' WDW equation, with 
the standard (Klcin-Gordon-like) invariant inner product 

($!, * 2 ) = % f dE a *J d a $ 2 (5) 

where the integral is to be taken over a spacelike hyper- 
surface inside the forward lightcone in /3-space. 

In order to construct solutions we separate variables 
by means of the ansatz $(p, z) = R(p)F(z) @, HH. For 
any eigenfunction F(z) obeying 

- A LB F(z) = EF(z) (6) 

the associated radial equation is solved by 

R ± (p)=p-^e ± ^ E -(^ 2 ^". (7) 

Positive frequency waves emanating from the singularity 
correspond to R-(p) and have positive inner product ([5]). 
To study the eigenvalues of the Laplace-Beltrami oper- 
ator on the unit hyperboloid we use a generalized upper 
half plane model z — (u, v) for the unit hyperboloid with 
co-ordinates u £ IR d ~ 2 and v £ R>o and the Poincare 
metric ds 2 = v~ 2 (dv 2 + du 2 ) such that 

A LB =v d - 1 d v (v 3 - d d v )+v 2 dl. (8) 

For the spectral problem we must specify boundary con- 
ditions. For the cosmological billiard, these are provided 
by infinite ('sharp') potential walls which encapsulate the 
effect of spatial inhomogeneities and matter fields near 
the spacelike singularity, as explained in [J, [5[ . Following 
the original suggestion of @, we are thus led to impose 
the vanishing of the wavefunction on the boundary of 
the fundamental domain specified by these walls. Ac- 
cordingly, let F{z) be any function on the hyperboloid 
satisfying ^ with Dirichlet conditions at the boundaries 
of this domain [2l[ . A direct generalization of the argu- 
ments on page 28 of Ref. [22j gives 

- (A LB F,F) > J dvd d - 2 uv 3 - d (d v F) 2 (9) 

with J6J and ([8|). Considering also 

(F,F) = J dvd^uv^F 2 

= -^-^ J dvd d - 2 uv 2 - d Fd v F , (10) 



the Cauchy-Schwarz inequality entails 

E >.(^)\ a.) 

From the explicit solution (J7J) we thus conclude that 
R±(p) — > when p — > 00, and therefore the full wave- 
function and all its p derivatives tend to zero near the 
singularity. While this conclusion also holds with Neu- 
mann boundary conditions (for which E > 0), the above 
inequality furthermore ensures that the full wavefunction 
is generically complex and oscillating. 

Let us now return to maximal supergravity, where the 
wavefunctions can be further analyzed by exploiting the 
underlying symmetry encoded in the Weyl group W(Eiq) 
and its arithmetic properties, and in particular the new 
links between hyperbolic We yl g roups and generalized 
modular groups uncovered in [231 ]. The Weyl reflections 
that the classical particle is subjected to when colliding 
with one of the walls are norm preserving, and there- 
fore the reflections can be projected to any hyperboloid 
of constant p, inducing a non-linear action on the co- 
ordinates z (given in (| 1 3|) below for the fundamental re- 
flections). For physical amplitudes to be invariant under 
the Weyl group, the full wavefunction (in /3-space) must 
transform as follows 

&(j>,z)=±*(fi,w I -z) (12) 

for the ten generating fundamental reflections Wj of 
W(E W ), labeled by I = -1,0, 1, ... ,8. Restricting the 
wavefunction to the fundamental Weyl chamber, one eas- 
ily checks that the plus sign in (fT2|) corresponds to Neu- 
mann boundary conditions, and the minus sign to Dirich- 
let conditions (which we adopt here), ^(p, z) is thus in- 
variant under even Weyl transformations s £ W + (Eiq) 
irrespective of the chosen boundary conditions. 

Choosing coordinates as in ([5]) the relevant variables 
now live in a 9-dimensional 'octonionic upper half plane' 
with z = u + iv where u = u £ O is an octonion. The ten 
fundamental reflections of W(Eiq) act as (j = 1, . . . , 8) 

w_i(z) = — , wq(z) = —6z9 + 6, uij(z) = —EjZSj (13) 

where z := u — iv, with iu = ui in accordance with 
Cayley-Dickson doubling [24[. Sj and 9, respectively, are 
the simple roots and the highest root of E s expressed 
as unit octonions [23|, [24| • We note that for u £ K the 
formulas (| 1 3[) reduce to the ones familiar from complex 
analysis, with z t— > 1/z, z 1— » — z + 1 and z ^ — z, 
generating the group PGL2(Z). For even Weyl trans- 
formations, we re-obtain the standard modular group 
PSL<i(7L) generated by S(z) — (w-iwi)(z) = —1/z and 
T(z) — (wqwi)(z) = z + 1. Similarly, the even Weyl 
group W + (Ei ) is isomorphic to PSX 2 (0) where are 
the integer octonions ('octavians', see (24|) and the group 
PSL 2 (Q) is defined by iterating the action of (|T5|) an 
even number of times [23| . Consequently, for maximal 
supergravity the bosonic wavefunctions ^(p, z) are odd 



3 



Maass wave forms for PS'Z^O), that is, invariant eigen- 
functions of the Laplace-Beltrami operator transforming 
with a minus sign in (| 12(1 under the extension W(E±q) 
of PSL2(0). Understanding these modular functions re- 
mains an outstanding mathematical challenge, see [25j 
for an introduction (and [13] for the PSL2CZ) theory). 
For the groups PSL2CZ) and PSL2(Z[i\) the (purely dis- 
crete) spectra of odd Maass wave forms have been inves- 
tigated numerically in [26|, [27], [H, [29|, [3(| ■ By modular 
invariance, the wavefunctions can be restricted to the 
fundamental domain of the action of W(E W ), and con- 
versely, their modular property defines them on the whole 
hyperboloid. The Klein-Gordon inner product ([5]) must 
likewise be restricted to the fundamental chamber 

($1, $ 2 ) = * J dvol(z)p d - 1 $t d p $ 2 . (14) 

where T is the intersection of C with the unit hyper- 
boloid. This is the canonical quantum gravity analog of 
the one-loop amplitude in string theory which is rendered 
finite upon 'division' by the modular group PSZ^Z). 

We now turn to the extension of the quantum billiard 
analysis to maximal supergravity and restrict to d = 10 
henceforth [3l], HH, [33| . Classically, the gravitino ?/v of 
D = 11 supergravity performs a separate fermionic bil- 
liard motion [34| . This is most easily expressed in a su- 
persymmetry gauge ipt — T t T a tp a (3l| and in the variables 

[13 (with r* = r 1 • • • r 10 ) 

<p a = g 1/4 T*T a iP a (no sum on a = 1,..., 10). (15) 

Using (|15p the Dirac brackets between two gravitino vari- 
ables (see (6.3) in [33|]) become y«} = —2iG ab 5 a [3, 
where we have written out the 32 real spinor components 
using the indices a, (3. The fermionic and bosonic vari- 
ables are linked by the supersymmetry constraint 

10 10 

S a EE ]T p a G aW b a = J2 *a<P% = . (16) 
a.b—1 a—1 

The supersymmetry constraint implies the Hamiltonian 
constraint H — by closure of the algebra 

\{S a ,Sp} = 8 aP H. (17) 

In order to quantize this system we rewrite the 320 real 
gravitino components in terms of 160 complex ones, 
and replace the Dirac brackets by canonical anticommu- 
tators to obtain a fermionic Fock space of dimension 2 160 
over the vacuum |0). For the supersymmetry constraint 

this amounts to the redefinition Sa = Sa + iSa+ig f° r 
A = 1, 16. The quantum constraint is then solved by 

16 

l*> = II & A ($( P ,z)\nj), (is) 

A=l 

(with <Sa|^) = 0) if and only if the function $(p, z) is 
a solution of the WDW equation 7Y$ = 0. While this 



solution is close to the 'bottom of the Dirac sea', there is 
an analogous one 'close to the top' with S\ replaced by 
Sa and \Q) by the completely filled state. 

The cosmological billiards description is very useful 
but takes into account the dependence on spatial inho- 
mogeneities and matter degrees of freedom only in a very 
rudimentary way via the infinite potential walls. It would 
thus be desirable to develop an approximation scheme for 
the quantum state in line with the 'small tension' expan- 
sion proposed in [ijj, and thereby hopefully resolve the 
difficulties encountered in extending the 'dictionary' of 
[l7T | to higher order spatial gradients and heights of roots 
in a quantum mechanical context. In the BKL approxi- 
mation, the full wavefunction is expected to factorize as 

i*mi> ~ n i**> ' ( i9 ) 

X 

near the singularity into a formal product over wave- 
functions of the type (Tj"5)) . one for each spatial point 
(with independent bosonic wavefunctions <& x (p(x), z(x)) 
and space- dependent diagonal metric variables /3 a (x) = 
(p(x), z(x)). The task is then to replace the formal ex- 
pression (|19p by a wavefunction depending on the (in- 
finite) tower of Eiq degrees of freedom, effectively im- 
plementing the de-emergence of space and time near 
the cosmological singularity, and their replacement by 
purely algebraic concepts [3l|, IHl- We note that as a 
consequence of the uniqueness of the standard bilinear 
(Cartan-Killing) form on Eiq there is a unique E±q ex- 
tension of the billiard Hamiltonian (0| given by 

mult (a) 

n ^ n + E E e- 2aW <s- (20) 

aeA + (Ei ) s=l 

where the first sum runs over the positive roots a of E\q . 
This extended system requires additional constraints. A 
first step in this direction was taken in [36| where a corre- 
spondence was established at low E\q levels between the 
classical canonical constraints of D = 11 supergravity (in 
particular, the diffeomorphism and Gauss constraints) on 
the one hand, and a set of constraints that can be consis- 
tently imposed on the E w / K{E W ) coset space dynamics 
on the other. The fact that the latter can be cast in 
a 'Sugawara-like' form as quadratic expressions in terms 
of the Eiq Noether charges [36] would make them par- 
ticularly amenable for the implementation on a quan- 
tum wavefunction. In addition, one would expect that 
PSL-2(Q) must be replaced by a much larger 'modular 
group' whose action extends beyond the Cartan subal- 
gebra degrees of freedom all the way into £10, perhaps 
along the lines suggested in (37j • 

As noted above, the inequality (jlip implies that 
$(/?, z) — > for p — > 00, and hence the wavefunction 
^ vanishes at the singularity, in such a way that the 
norm is preserved. Its oscillatory nature entails that it 
cannot be analytically extended beyond the singularity, 
a result whose implications for the question of singular- 
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ity resolution in quantum cosmology remain to be ex- 
plored. The mechanism usually invoked to resolve sin- 
gularities in canonical approaches to quantum geometro- 
dynamics would be to replace the classical 'trajectory' in 
the moduli space of 3-geometries (that is, WDW super- 
space) by a quantum mechanical wave functional which 
'smears' over the singular 3-geometries. By contrast, the 
present work suggests a very different picture, namely the 
'resolution' of the singularity via the effective disappear- 
ance (de-emergence) of space-time near the singularity 
(see also [HI). The singularity would thus become effec- 
tively 'unreachable'. This behavior is very different from 
other possible mechanisms, such as the Hartlc-Hawking 
no boundary proposal [38[, or cosmic bounce scenarios 
of the type considered recently in the context of minisu- 
perspace loop quantum cosmology [H, Sl| , both of 
which require continuing the cosmic wavepacket into and 
beyond the singularity at p = oo. 

A key question for singularity resolution concerns the 



role of observables, and their behavior near the singu- 
larity While no observables (in the sense of Dirac) are 
known for canonical gravity, we here only remark that for 
the Eio/K(Eio) coset model the conserved E\o Noether 
charges do constitute an infinite set of observables, as 
these charges can be shown to commute with the full 
Eio Hamiltonian (f!?0|) . The expectation values of these 
charges are thus the only quantities that remain well- 
defined and can be sensibly computed in the deep quan- 
tum regime, where the E w / K(E W ) coset model is ex- 
pected to replace space-time based quantum field theory. 
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